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Abstract 

As  an  example  of  a singular  singular- 
perturbation  problem,  we  solve  the  high  gain 
systea  x * A(t)x  + B(t)u,  u » gC(t)x  as 
g - 1/u  ■»  Our  primary  assumption  requires 
BC  to  have  fixed  rank  k < n and  no  unstable 
eigenvalues.  We  find  an  asymptotic  solution  of 
the  form  x(t,u)  * X(t,u)  + H(t,u)  where  the 
fast  transient  E ■»  0 as  T * t/w  -*  ". 


We  shall  consider  the  linear  high  gain  system 

Ix  • A(t)x  + B(t)u,  t ^ 0 
u - gC(t)x 


(1) 


where  the  state  vector  x is  n-dinensional,  the 
control  vector  u is  m-dinensional , and  the  sca- 
lar gain  factor  g is  large.  Introducing  the 
small  parameter 


I/g. 


(2) 


- 1.  Introduction 


then,  the  feedback  system  (1)  takes  the  form 


Singular  perturbation  methods  have  now  become 
an  increasingly  accented  and  valued  new  approach 
for  solving  limiting  problems  an-  obtaining  sim- 
plified models  in.  central  ic Kokotovic  et  al. 
(1976)).  Recent!. •,  Jameson  ur.d  this  author  [cf. 
O'Malley  and  Jameses  (1975,  1 -•  7 T ) j have  .-lined  new 
insight  on  singular  problems  cf  :pt;-al  control  by 
analyzing  them  as  the  lima:  of  r.-arly  singular 
ones.  This  involved  the  saucy  of  singularly  per- 
turbed tvo-geini  boundary  val orci'.vr.s.  Unlike 
the  more  traa.tlor.al  trchlems  (ii  =-.cauntcrcd  with 
linear  regulator  problems  with  singularly  perturbed 
state  equations),  the  limiting  be-.avior  is  not  de- 
termined by  the  solution  of  a natural  reduced  prob- 
lem obtained  by  setting  c » 3.  Instead,  the  re- 
duced problem  c.as  an  infinite  number  of  solutions 
and  some  ocher  means  is  r.eedec  to  specify  which  of 
these  is  the  limiting  solution.  Motivated  by  this 
problem,  and  the  critical  need  to  sclve  analogous 
problems  for  stiff  differential  equations,  O'Malley 
and  Flaherty  (1977)  have  begun  to  examine  such 
"singular”  singular-percurfcation  problems.  They 
turn  out  to  be  interesting  mathematically,  in  addi- 
tion to  their  obvious  practical  importance. 


In  these  proceedings,  Young  et  ai.  (1977)  dis- 
cuss high  gain  feedback  systems.  We  refer  the 
reader  to  their  paper  for  appropriate  references 
and  a discussion  of  the  significance  of  such  stud- 
ies in  various  control  situations.  We  shall  exam- 
ine these  as  an  example  of  a singular  sinqular- 
,perturbation  problem  illustrating  a very  successful 
interplay  between  singular  perturbations  theory  and 
its  applications  in  control. 
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MX  - (B(t)C(t)  + pA(t))x,  t>0.  (3) 

If  BC  remains  a stable  matrix,  classical  singu- 
lar perturbations  theory  implies  the  existence  of 
a unique  asymptotic  solution  to  the  initial  value 
problem  for  (3)  of  the  form 


x(t,u) 


» n(i,u) 


"V  E 

J-0 


iyou3. 


(*> 


Hero,  the  terms  H 
stretched  variable^ 


each  tend  to  zero  as  the 


T - t/p 


(S) 


tends  to  infinity  (cf.  Hoppensteadt  (1971)  and 
O'Malley  (1974)).  In  particular,  away  from  t * 0, 
the  limiting  solution  X (t)  * 0 Is  the  unique 
solution  of  the  reduced  ° problem 


B(t)C(t)X  (t)  =>  0.  (6) 

o 

Young  et  al.  consider  the  singular  problem  where 
rank  (BC)  * ci  < n.  Then,  (6)  has  an  n - n dim- 
ensional linear  manifold  of  solutions,  though  ve 
expect  the  initial  value  problem  to  have  a unique 
solution  for  each  u > 0. 


We  shall  make  the  assumption 

{the  n x n matrix  B(t)C(t)  has  a con- 
stant rank  k,  0 £ k < n,  for  all  t 
In  0 ^ t £ T < »;  its  null  space  is 
spanned  by  n - k linearly  independent 
eigenvectors;  and  its  nonzero  eigenval- 
ues all  have  negative  real  parts  there. 
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Then,  we  shall  obtain  a unique  asymptotic  solution 
to  the  initial  value  problem  for  (3)  of  the  form 


guaranteed  by  hypothesis  (H) 


x(t,p)  - X(t,u)  + n(r,p) 


3.  A Solution  via  a Transformed  Problem 


where  the  outer  solution  X(C,u)  and  the  boundary 
layer  correction  E(t,p)  both  have  asymptotic 
series  expansions  in  u,  with  the  boundary  layer 
terms  decaying  to  zero  as  the  stretched  variable 
T » t/p  tends  to  infinity. 


If  we  make  the  1-1  transformation 

w » Ex,  (1 

the  variable  w satisfies  the  singularly  per- 
turbed system 


2.  Preliminary  Linear  Algebra 


UV  - (EBCE’  p(E  + EA)E')w. 


Under  hypothesis  (H),  the  matrix  BC  can  be 
put  into  Its  row-reduced  echelon  form  by  use  of  an 
orthogonal  macrix  E,  i.e.  we'll  obtain 


This,  however.  Is  also  a singular  singular-per- 
turbation problem  since  the  matrix  E3CE'  has 
rank  k < n. 


E(t)B(t)C(t) 


where  U(t)  is  a k x n matrix  of  rank  k.  In- 
deed, E can  be  readily  obtained  via  Householder 
transformations  (cf.  Colub  (1965)).  Writing 


where  E^  Is  k x n,  we  have 
E2EC  - 0 

while  the  orthogonality  of  E implies  that 


The  structure  of  the  solution  becomes  obvi- 
ous upon  splitting  w into  vector  components 


They  satisfy  the  "nonsingular"  singularly  per- 
turbed system 


Ww^  - Sw^  + EjBCE^  + + pfr^w. 

w m E w + F v 
2 3 1 T2 


where  we've  used  the  decomposition  (14)  for  EBCE' 
and  set 


E F'  « 0 rr’«I 
fcl  2 ’ “1*1  k’ 

Further, 


-vi  *Vk-  {11> 


? - E^E,  - ? and  Q - E^  - Q*  (1) 

are  complementary  projections  with  ? + Q * I_, 

and  Q projects  into  the  roil  space  of  (BC)'" 
since 

(BC) 'Q  - (3C)'E^E2  - (E;3C;'E2  = 0 

by  (10).  Since  the  dimensions  of  R(Q),  the 
range  of  Q,  and  M((EC)'),  the  null  space  of 
(BC) ' , are  both  n - k,  we  have 


R(Q)  " •'•'((BC)’) 


R(P)  - R(BC). 


Finally,  (9)  and  (10)  imply  that 


EjBCE^  EjBCE^ 


so  the  k x k matrix 

S - E^CE’  is  stable  ( 

since  It  has  the  k stable  eigenvalues  of  BC 


(E  + EA)E‘ 


Standard  singular  perturbation  theory  shows  that 
the  initial  value  problem  for  (19)  will  have  a 
unique  solution  of  the  form 

jw1(t,y)  - wx(t,u)  + Pj^TpU) 


Jv^Ct.u)  * W2(e,p)  + wP2(r,ii) 

where  all  term3  have  asymptotic  expansions  in  u 
and  P and  P?  tend  to  zero  as  T •»  <».  Fur- 
ther, the  leading  terms  and  W of  the 

outer  expansion  satisfy  the  reduced  problem. 

( 0 - SM10  + ElB<iEiW20 


W20  * E3W10  + E4M20’  W20(0)  " U2(0)  ” E2(0)x(0)' 


W10  * -S'lElBCEiM20 

while  W2q  uniquely  satisfies 

*20  ’ (E4  - E3S'1BCE’)W20. 

W20(0)  ■ V0)*(°>- 


■ life:  TTT 


The  leading  terms  of  the  boundary  layer  correction 
satisfy 


dP 


10 


dt 


s<0)p10.  " E,(0)x(0)  - W,n(0) 


10 


10' 


^ • EJ<">V 


Thus,  the  decaying  solutions  are  given  by 
fpi0(T)  - eS(0)T(Ei(0)x(0)  - W10(0)) 


P20^  -fx 


C3(0)P10(s)ds 


E3(0)S-1(0)P10(t). 


(25) 


Higher  order  terms  follow  analogously.  Finally, 

’ v 

X * E’  { | implies  that  the  unique  solution  of 


the  original  Initial  value  problem  Is  represented 
In  the  additive  fora  (7)  with  the  outer  solution 

X(t,u)  - E’Ct'.ulV^t.u)  + E^(t,u)W2(t.u)  (26) 

and  the  boundary  layer  correction 

n(T,u)  - E*(ct, (r,u) 

+ uE^(£t.u)?2(t,--) 


(27) 


where  B is  the  projection  B = I - AeC.  Thus 
each  X has  been  successively  determined  up  to 
Its  projection  QX^  onto  N((BC)'). 

Since  QBC  - 0,  consistency  of  (30)  requires 
that  QE  3 - 0.  Since  E ” 0,  this  Fredhola 
alternative  argument  requires  that  (QX  )" 


» (Q  + QA)X  for  each  j _>  0.  Using 


y 


(32), 


V X • O ' * 

then,  we  finally  obtain  the  linear  differential 
equation 


(QXj)' 


(Q  + QA)B(QX^)  + (Q  + QA)AEj_1 


(33) 


for  QX.  • Ve  therefore  find  the  outer  expansion 
by  solving  algebraic  equations  (31)  for  PX  and 
the  differential  equation  (33)  for  Its  complement 
QX  . Thus,  we're  able  to  formally  obtain  X(c,u) 
termwise  up  to  specifying  the  initial  values 
Q(0)X  (0).  Since  the  initial  value  x(0)  won't 
generally  satisfy  (31)  at  t - 0,  we'll  need  the 
boundary  layer  correction  H to  account  for  the 
resulting  initial  jump  (l.e.  a nonuniform  conver- 
gence at  t » 0 occurring  in  R(BC)). 


By  linearity,  n(t,u)  must  satisfy 
(B(et)C(ct)  + uA(et)]n 
as  a power  series 


dll 

dt 


ji(t,p)  “v.  z n.(t)u 


i 


J-0 


r 


(34) 


(35) 


In  p.  Thus,  we  must  successively  obtain  decaying 
solutions  of 


which  decays  to  zero  as  t — 


4.  A Direct  Solutitr. 

Let  us  now  try  to  directly  obtain  a solution 
of  the  form  (7).  Since  tne  bcur.C3ry  layer  correc- 
tion ,1  decays  to  zero  as  - - »,  it's  necessary 
that  the  outer  solution  X satisfy  the  system  (3) 
for  all  C > 0.  Thus  the  cuter  expansion 


x(t,u)  n.  r r (t)J 

J-o  j 


(28) 


must  satisfy 

lii  - (BC  + uA)X  (29) 

as  a power  series  in  p . This  requires  that 

®CXJ  ‘ Xj-1  ' ^j-l  5 5j-l  (30) 

for  each  j > 0 with  E_^  S 0.  Multiplying  this 
equation  by  E^  and  manipulating,  we  obtain 

PXj  - -AbC(QXj)  + Ae^  t (31) 

for  A » E'S  ^E, . But  P + Q * I then  implies 
that  11  n 


XJ  “ + A5J-1 


(32) 


dn 

tf-  - B(o)c(0)n  + c , 


i > o. 


(36) 


where  E ^ = 0 and  E,  , is  generally  a linear 
combination  of  preceding  terms  IT  with  polynomi- 
al coefficients  in  T.  Since  QBC  * 0,  it  fol- 


lows that  — (Q(0)TT  (t ) ) 
dt  J 


Q(0)Ej_i(t),  so 


Q(0)nj(T)  - -Q(0)  j cj_1(s)ds,  J > 0 


(37) 


It  remains  to  find  the  complementary  P(0)IT  (t)'s. 
Multiplying  (36)  by  E^(0),  we  have  ' 


<Ei(o>y 


S(0)(E1(0)Bj) 


+ E1(0)(B(0)C(0)Q(0)Hj(t)  + C x). 


integrating,  then,  and  using  (37),  we  obtain 
lyr)  - E^(0)eS(0)TE1  (0)11^0) 

j E^_1(s)ds 

qco,  f 


- Q(0) 


+ E 


es(°)(T-r)Ei(0). 


. i _J 


•(-B(0)C(0)Q(0)  C (s)ds 
■'r  J 

+ CJ_1(r))dr.  (38) 

Since  C = 0,  II0(t)  is  determined  up  to  its 
initial  value  E ^ ( 0) (0) * and  it  is  exponential- 
ly decaying.  Using  induction,  it  readily  follows 
that  succeeding  terms  H (t ) are  exponentially 
decaying  whatever  initial  values  are  later  selec- 
ted. 

To  completely  specify  the  expansions  genera- 
ted, we  must  obtain  the  initial  values  Q(Q)X.(0) 
and  E (0)1!  (0)  termwise.  Since  (37)  implies 
that  Q(0)nJ (0)  * 0,  the  representation  (7)  1m- 
plies  that  °Q(0) X (0)  = Q(0)x(0)  and  the  limit- 
ing outer  solution  X will  be  completely  deter- 
mined a3  the  solution  of  the  n - k th  order 
dynamical  system 


mlc  system  (39)  In  an  n - k dimensional  space 
where  k ■ rank  (BC).  The  results  continue  to 
hold  for  all  t 0 provided  (H)  remains  valid 
and  X (t)  decays  exponentially  to  zero  as 
t ■*  “,°  Finally,  observe  that  (44)  displays  the 
two  time  scale  nature  of  the  limiting  solution  on 
0 < t < T. 

Our  results  should  be  further  related  to 
earlier  work  in  the  literature.  We  note,  in  par- 
ticular, Young  et  al.’s  result  that  CX  * 0 
when  k ■ m “ rank  (CB).  It  Implies  that  the 
control  would  stay  bounded  away  from  t * 0. 


Acknowledgment:  Thanks  should  be  given  to  Petar 

Kokocovic  and  his  co-authors  for  emphasizing  the 
singular  perturbation  nature  of  large  gain  sys- 
tems. It  fit  well  with  the  Cheory  of  singular 
singular-perturbation  problems,  currently  under 
study  by  this  author. 


(QX  r » (Q  + QA)B(QX  ), 
o o 

Q(0)X  (0)  - Q(0)x(0) 

O 


(39) 


together  with 

• X - B(QX  ).  (40) 

o o 

Then,  £.(0)1!  (0)  - E_(0)x(0)  - E.(0)X  (0)  and 
lo  1 1 o 

(37)  implies  tha; 

n (T)  - E’(C)eS'0)T  E.(0)(x(0)  - X (0)).  (41) 

O 1 1 o 


Continuing  for  each  j > 0,  we  sec 

3(C)X  (0)  - -Q(0)r,(0) 

- Q(0)  ; i.  . (s) d s 
J 3 J_1 


(42) 


and  this  allows  us  to  uniquely 
wise.  Then 


obtain  X^ 


term- 


Ex (0)^(0)  - -E1(0)X.(0) 


(43) 


becomes  known  and  we  completely  find  the  boundary 
layer  correction  term  H . The  solution  generated 
will,  of  course,  agree  ^ with  that  obtained  via 
the  preceding  transformed  problem. 


5.  Summary 


We  find  that  the  asymptotic  solution  of  the 
high  gain  feedback  problem  (1),  i.e.,  the  Initial 
value  problem  for  (3)  is  of  the  form 

x(t,u)  - XQ(t)  + IIo(t/p)  + 0(u)  (44) 

under  hypothesis  (H) . Specifically,  we  note  that 
the  fast  transient  It  decays  to  zero,  remaining 
in  R(BC),  while  the°  outer  limit  X satisfies 
BCXq  ■ 0.  Moreover,  Xq  is  determined  by  a dyna- 
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